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ABSTRACT 

Context. In the Grand Tack (GT) scenario for the young solar system, Jupiter formed beyond 3.5 AU from the Sun and migrated 
as close as 1.5 AU until it encountered an orbital resonance with Saturn. Both planets then supposedly migrated outward for several 
10^ yr, with Jupiter ending up at» 5 AU. The initial conditions of the GT and the timing between Jupiter’s migration and the formation 
of the Galilean satellites remain unexplored. 

Aims. We study the formation of Ganymede and Callisto, both of which consist of50% H 2 O and rock, in the GT scenario. We 
examine why they lack dense atmospheres, while Titan is surrounded by a thick N 2 envelope. 

Methods. We model an axially symmetric circumplanetary disk (CPD) in hydrostatic equilibrium around Jupiter. The CPD is warmed 
by viscous heating, Jupiter’s luminosity, accretional heating, and the Sun. The position of the H 2 O ice line in the CPD, which is 
crucial for the formation of massive moons, is computed at various solar distances. We assess the loss of Galilean atmospheres due to 
high-energy radiation from the young Sun. 

Results. Ganymede and Callisto cannot have accreted their H 2 O during Jupiter’s supposed GT, because its CPD (if still active) was 
too warm to host ices and much smaller than Ganymede’s contemporary orbit. From a thermal perspective, the Galilean moons might 
have had significant atmospheres, but these would probably have been eroded during the GT in < 10^ yr by solar XUV radiation. 
Conclusions. Jupiter and the Galilean moons formed beyond 4.5 ± 0.5 AU and prior to the proposed GT. Thereafter, Jupiter’s 
CPD would have been dry, and delayed accretion of planetesimals should have created water-rich lo and Europa. While Galilean 
atmospheres would have been lost during the GT, Titan would have formed after Saturn’s own tack, because Saturn still accreted 
substantially for » 10® yr after its closest solar approach, ending up at about 7 AU. 

Key words. Accretion, accretion disks - Planets and satellites: atmospheres - Planets and satellites: formation - Planets and satellites: 
physical evolution - Sun: UV radiation 


1. Introduction 


Recent simulations of the early solar system suggest that the four 
giant planets underwent at least two epochs of rapid orbital evo¬ 
lution. In the Grand Tack (GT) model (Walsh et al. 2011), they 
formed beyond the solar water (H 2 O) ice line (Ciesla & Cuzzi 
2006) between 3.5 AU and 8 AU. During the first few 10® yr in 
the protostellar disk, when Jupiter had fully formed and Saturn 
was still growing, Jupiter migrated as close as 1.5 AU to the Sun 
within some 10® yr, until Saturn gained enough mass to migrate 
even more rapidly, catching Jupiter in a 3:2 or a 2:1 mean motion 
resonance (Pierens et al. 2014). Both planets then reversed their 
migration, with Jupiter ending up at approximately 5 AU. The 
other gas planets also underwent rapid orbital evolution due to 
gravitational interaction, pushing Uranus and Neptune to and be¬ 
yond about 10 AU, respectively, while Saturn settled at » 7 AU. 
A second period of rapid orbital evolution occurred « 7 x 10*^ yr 
later, according to the Nice model (Tsiganis et al. 2005; Mor- 
bidelli et al. 2005; Gomes et al. 2005), when Jupiter and Saturn 
crossed a 2:1 mean motion resonance, thereby rearranging the 
architecture of the gas giants and of the minor bodies. 
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While the GT delivers adequate initial conditions for the 
Nice model, the initial conditions of the GT itself are not well 
constrained (Raymond & Morbidelli 2014). Details of the mi¬ 
grations of Jupiter and Saturn depend on details of the solar ac¬ 
cretion disk and planetary accretion, which are also poorly con¬ 
strained (Jacobson & Morbidelli 2014). We note, however, that 
the moons of the giant planets provide additional constraints on 
the GT that have scarcely been explored. As an example, using 
A-body simulations, Deienno et al. (2011) studied the orbital 
stability of the Uranian satellites during the Nice instability. Hy¬ 
pothetical moons beyond the outermost regular satellite Oberon 
typically got ejected, while the inner moons including Oberon 
remained bound to Uranus. Their results thus support the valid¬ 
ity of the Nice model. Later, Deienno et al. (2014) found that 
one of three Jovian migration paths proposed earlier (Nesvorny 
& Morbidelli 2012) is incompatible with the orbital stability and 
alignment of the Galilean moons, yielding new constraints on 
the Nice model. 

To our knowledge, no study has used the Galilean moons or 
Saturn’s major moon Titan, which is surrounded by a thick ni¬ 
trogen (N 2 ) atmosphere, to test the plausibility of the GT model. 
In this Letter, we focus on the early history of the icy Galilean 
moons, Ganymede and Callisto, and on Titan in the GT scenario; 
assuming that the GT scenario actually took place in one form or 


Article number, page 1 of 4 




A&A proofs: manuscript no. ms 


another, we identify new constraints on the timing of their for¬ 
mation. The novel aspect of our study is the evolution of the H 2 O 
ice line in Jupiter’s circumplanetary disk (CPD) under the effect 
of changing solar illumination during the GT. As the position 
and evolution of the ice line depends on irradiation from both 
the forming Jupiter and the Sun (Heller & Pudritz 2015a), we 
want to understand how the observed dichotomy of two mostly 
rocky and two very icy Galilean moons can be produced in a GT 
setting. 

2. The icy Galilean moons in the Grand Tack model 

2.1. Formation in the circumjovian accretion disk 

Both Ganymede and Callisto consist of about 50% rock and 
50 % H 2 O, while the inner Galilean satellites lo and Europa are 
mostly rocky. This has been considered a record of the tempera¬ 
ture distribution in Jupiter’s CPD at the time these moons formed 
(Pollack & Reynolds 1974; Sasaki et al. 2010). In particular, the 
H 2 O ice line, which is the radial distance at which the disk is 
cool enough for the transition of H 2 O vapor into solid ice, should 
have been between the orbits of rocky Europa at about 9.7 Jupiter 
radii (Rjup) and icy Ganymede at about 15.5 Rjup from Jupiter. 

We simulate a 2D axisymmetric CPD in hydrostatic equilib¬ 
rium around Jupiter using a “gas-starved” standard disk model 
(Canup & Ward 2002, 2006) that has been modified to include 
various heat sources (Makalkin & Dorofeeva 2014), namely, (i) 
viscous heating, (ii) planetary illumination in the “cold-start sce¬ 
nario”, (iii) direct accretion onto the CPD, and (iv) stellar illu¬ 
mination. The model is coupled to pre-computed planet evolu¬ 
tion tracks (Mordasini 2013), and we here assume a solar lu¬ 
minosity 0.7 times its current value to take into account the 
faint young Sun (Sagan & Mullen 1972). Details of our semi- 
analytical model are described in Heller & Pudritz (2015a,b). 

We evaluate the radial position of the circumjovian H 2 O ice 
line at different solar distances of Jupiter and study several disk 
opacities (/cp) as well as different shutdown accretion rates for 
moon formation (Mshut), all of which determine the radial dis¬ 
tributions of the gas surface density and midplane temperature 
in the CPD. We follow the CPD evolution after the planet opens 
up a gap in the circumstellar disk (CSD) until the planetary ac¬ 
cretion rate (Mp) in the pre-computed track drops to a particu¬ 
lar value of Mshut- In Heller & Pudritz (2015a,b) we found that 
Afshut = lOMcanMyr^* positions the H 2 O ice line between Eu¬ 
ropa and Ganymede for a broad range of Kp values. 

To assess the effect of heating on the position of the H 2 O ice 
line in the CPD as a function of solar distance, we compare two 
CSD models. Eirst, we use the model of Hayashi (1981, H81, his 
Eq. 2.3), which assumes that the CSD is mostly transparent in the 
optical. Second, we use the model for an optically thick disk of 
Bitsch et al. (2015, B15, their Eqs. A.3 and A.7), which takes 
into account viscous and stellar heating as well as radiative cool¬ 
ing and opacity transitions. We consider both a low-metallicity 
(Z = 0.001) and a solar-metallicity star (Z = 0.02) with a stel¬ 
lar accretion rate M* = 3.5 x 10“®Moyr“^ corresponding to 
Jupiter’s runaway accretion phase (Mordasini 2013), in which 
the planet opens up a gap and moon formation shuts down. 

Eigure 1 shows the results of our calculations, assuming fidu¬ 
cial values/cp = 10“^m“^kg“* and Mshut = lOMcanMyr^* as 
well as a centrifugal CPD radius as per Machida et al. (2008). 
The abscissa denotes distance from the young Sun, the ordinate 
indicates distance from Jupiter. The vertical dotted line high¬ 
lights the critical solar distance (flcrit) at which the CPD loses 
its H 2 O ice line in the Z = 2% H15 model. At 5.2 AU, Cal- 
listo’s position (circle labeled “C”) at the outer CPD edge is due 
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Fig. 1. Radial distance of Jupiter’s H 2 O ice line (roughly horizontal 
lines) and centrifugal disk radius (curved solid line) during its final ac¬ 
cretion phase. Inside about 4.5 AU from the Sun, Jupiter’s disk does not 
contain an H 2 O ice line in any of the CSD models (Hayashi 1981; Bitsch 
et al. 2015, see legend). Locations of the Galilean moons are indicated 
by symbols at 5.2 AU. Symbol sizes scale with the physical radii of the 
moons (orange; rocky, blue: icy composition). 

to our specific CPD scaling (Heller & Pudritz 2015b), whereas 
Ganymede’s position (circle labeled “G”) near the H 2 O ice line 
is not a fit but a result.* The CPD radius shrinks substantially to¬ 
wards the Sun owing to the increasing solar gravitational force, 
while the ice line recedes from the planet as a result of enhanced 
stellar heating. A comparison of the different CSD models (see 
legend) indicates that variations of stellar metallicity or solar il¬ 
lumination have significant effects on the circumjovian ice line, 
but the critical effect is the small CPD radius in the solar vicinity. 

Most importantly. Pig. 1 suggests that Ganymede and Cal¬ 
listo cannot have accreted their icy components as long as Jupiter 
was closer than about aoit = 4.5 AU to the faint Sun, where the 
H 2 O ice line around Jupiter vanishes for both the H81 and the 
B15 model. We varied Kp and Mshut by an order of magnitude 
(not shown), which resulted in changes of this critical solar dis¬ 
tance of < 0.5 AU. Moreover, closer than 5.2 AU from the Sun, 
Callisto’s contemporary orbit around Jupiter would have been 
beyond the CPD radius, though still within Jupiter’s Hill sphere. 

Hence, large parts of both Ganymede and Callisto cannot 
have formed during the GT, which supposedly brought Jupiter 
much closer to the Sun than 4.5 + 0.5 AU. Thereafter, Jupiter’s 
accretion disk (if still present at that time) would have been 
void of H 2 O because water first would have vaporized and then 
been photodissociated into hydrogen and oxygen. If Ganymede 
or Callisto had acquired their H 2 O from newly accreted plan- 
etesimals after the GT (e.g. through gas drag within the CDP, 
Mosqueira & Estrada 2003), then lo (at 0.008 Hill radii, Rh) 
and Europa (at 0.01 Rh) would be water-rich, too, because plan- 
etesimal capture would have been efficient between 0.005 and 
0.01 Rh (Tanigawa et al. 2014). Hence, Ganymede and Callisto 
must have formed prior to the GT, at least to a large extent. 

2.2. Atmospheric escape from the Gaiiiean moons 

Not only must a successful version of the GT model explain the 
observed H 2 O ice contents in the Galilean moons, it also must be 


’ Heller & Pudritz (2015b) argue that this suggests that Jupiter’s H 2 O 
ice line acted as a moon migration trap for Ganymede. 
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Fig. 2. Integrated Maxwell-Boltzmann velocity distribution of N 2 and 
H 2 O gas molecules with temperatures similar to Callisto’s surface tem¬ 
perature during accretion. The black vertical line denotes Callisto’s es¬ 
cape velocity, values for the other Galilean moons are indicated by gray 
vertical lines. The black circle refers to an example discussed in the text. 

compatible with both the observed absence of thick atmospheres 
on the Galilean moons and the presence of a massive N 2 atmo¬ 
sphere around Saturn’s moon Titan. 

Several effects can drive atmospheric escape: (1) thermal (or 
“Jeans”) escape, (2) direct absorption of high-energy (X-ray and 
ultraviolet, XUV) photons in the upper atmosphere, (3) direct 
absorption of high-energy particles from the solar wind, (4) im¬ 
pacts of large objects, and (5) drag of heavier gaseous compo¬ 
nents (such as carbon, oxygen, or nitrogen) by escaping lighter 
constituents (such as hydrogen) (Hunten et al. 1987; Pierrehum- 
bert 2010).^ Although additional chemical and weather-related 
effects can erode certain molecular species (Atreya et al. 2006), 
we identify the XUV-driven non-thermal escape during the GT 
as a novel picture that explains both the absence and presence of 
Galilean and Titanian atmospheres, respectively. 

2.2.1. Thermal escape 

As an example, we assess Callisto’s potential to lose an initial 
N 2 or H 2 O vapor atmosphere via thermal escape. Our choice is 
motivated by Callisto’s relatively low surface temperature during 
accretion, which was about 300 K at most, taking into account 
heating within Jupiter’s accretion disk and from the accretion 
of planetesimals (Lunine & Stevenson 1982). The other moons 
might have been subject to substantial illumination by the young 
Jupiter, so the effect of thermal escape might be harder to assess. 

r-OO 

We compute Fme)*)) = J- dvfMB(v), with /mb)^) as the 
Maxwell-Boltzmann velocity distribution, for various gas tem¬ 
peratures V. Then T’^™(£)=T'MB(t')/7^MB(0) is an approximation 
for the fraction of the atmosphere that is above a velocity D. 

In Figure 2 we compare 7 ^m™( 0) for N 2 (black lines) and 
H 2 O (red lines) molecules at temperatures between 200 K and 
400 K (see legend) with Callisto’s gravitational escape veloc¬ 
ity (2.44 km s“*, black vertical line). As an example, T’M™(y > 
2.44 km s“') < 10“'"^ for an N 2 troposphere at 300 K (see black 
circle), suggesting that a negligible fraction of the atmosphere 
would be beyond escape velocity. 


^ Atmospheric drag (5) can only occur if any of the other effects (1)- 
(4) is efficient for a lighter gas component. 


Taking into account collisions of gas particles, their finite 
free path lengths, and the extent and temperature of the exo¬ 
sphere, Pierrehumbert (2010) estimates the loss time of N 2 from 
Titan as 10^"^ yr. For Callisto, which has about 80% of Ti¬ 
tan’s mass, this timescale would only be reduced by a factor of a 
few (Pierrehumbert 2010, Eq. 8.37) even if it had been as close 
as 1.5 AU from the Sun. At that distance, Callisto would have 
received about 40 times higher irradiation than Titan receives 
today (assuming an optically thin CSD), and its exosphere tem¬ 
peratures might have been about 40'^"^ 2.5 times higher than 

Titan’s. 

Hence, thermal escape cannot be the reason for the lack of 
an N 2 atmosphere on Callisto, even during the GT as close as 
1.5 AU from the Sun. 

2.2.2. XUV-driven non-thermal escape and drag 

The X-ray and UV luminosities of the young Sun were as high 
as 10^ (Ribas et al. 2005) and 10^ (Zahnle & Walker 1982) times 
their current values, respectively, raising the question whether 
direct absorption of high-energy photons or atmospheric drag 
would have acted as efficient removal processes of an early Cal- 
listonian atmosphere. Such a non-thermal escape would have 
eroded a hypothetical initial N 2 atmosphere from Earth in only a 
few 10® yr (Lichtenegger et al. 2010), owing to the high exobase 
temperatures (7000-8000K) and the significant expansion of 
the thermosphere above the magnetopause (Tian et al. 2008). 
However, the Earth’s primitive atmosphere was likely C02-rich, 
which cooled the thermosphere and limited the N 2 outflow. 

If Callisto initially had a substantial CO 2 or H 2 O steam at¬ 
mosphere, perhaps provided by outgassing, both gases would 
have been photo-dissociated in the upper atmosphere, which 
then would have been dominated by escaping H atoms; N 2 and 
other gases would have been dragged beyond the outer atmo¬ 
sphere and lost from the moon forever. Lammeret al. (2014) sim¬ 
ulated this gas drag for exomoons at 1 AU from young Sun-like 
stars and found that the H, O, and C inventories in the initially 
thick CO 2 and H 2 O atmosphere around moons 10% the mass 
of the Earth (four times Ganymede’s mass) would be lost within 
a few 10® yr depending on the initial conditions and details of 
the XUV irradiation. Increasing the moon mass by a factor of 
five in their computations, non-thermal escape times increased 
by a factor of several tens. Given that strong a dependence of 
XUV-driven escape on a moon’s mass, Ganymede and Callisto 
would most certainly have lost initial N 2 atmospheres during the 
GT (perhaps within 10^ yr), even if they approached the Sun as 
close as 1.5 AU rather than 1 AU as in Lammer et al. (2014). 

2.3. Titan’s atmosphere in the Grand Tack modei 

If XUV-driven atmospheric loss from the Galilean moons indeed 
occurred during the GT, this raises the question why Titan is 
still surrounded by a thick N 2 envelope, since Saturn supposedly 
migrated as close as 2 AU to the young Sun (Walsh et al. 201 1). 
We propose that the key lies in the different formation timescales 
of the Galilean moons and Titan. 

In the GT simulations of Walsh et al. (201 1), Saturn accretes 
about 10% of its final mass over several 10® yr after its tack, 
when Jupiter is already fully formed. Hence, while the Galilean 
moons must have formed before Jupiter’s GT (Sect. 2.1) and 
migrated towards the Sun, thereby losing any primordial at¬ 
mospheres, Titan formed after Saturn’s tack and on a longer 
timescale. Titan actually must have formed several 10® yr after 
Saturn’s tack, or it would have plunged into Saturn via its own 
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type I migration within the massive CPD (Canup & Ward 2006; 
Sasaki et al. 2010) because there would be no ice line to trap it 
around Saturn (Heller & Pudritz 2015b). 

The absence of N 2 atmospheres around the Galilean satellites 
supports the GT scenario and is compatible with the presence 
of a thick N 2 atmosphere around Titan. The Galilean satellites 
likely lost their primordial envelopes while approaching the Sun 
as close as 1.5 AU, whereas Titan formed after Saturn’s tack at 
about 7 AU under less energetic XUV conditions. Its N 2 atmo¬ 
sphere then built up through outgassing of NH 3 accreted from 
the protosolar nebula (Mandt et al. 2014). 

3. Discussion 

Although the GT paradigm is still controversial, our results pro¬ 
vide additional support for it. We show that the outcome of 
atmosphere-free icy Galilean satellites and of a thick N 2 atmo¬ 
sphere around Titan is possible in the GT scenario. This demon¬ 
strates how moon formation can be used to constrain the mi¬ 
gration and accretion history of planets - in the solar system 
and beyond. The detection of massive exomoons around the ob¬ 
served exo-Jupiters at 1 AU from Sun-like stars would indicate 
that those gas giants migrated from beyond about 4.5 AU. 

We tested two CSD models (H81 and B15) to study the 
effect of changing solar heating on the radial position of the 
H 2 O ice line in Jupiter’s CPD. Both models yield similar con¬ 
straints on the critical solar distance (c/cnt) beyond which the 
icy Galilean satellites must have formed. We varied opacities, 
shutdown accretion rates, and stellar metallicities to estimate the 
dependence of our results on the uncertain properties of the pro¬ 
toplanetary and protosatellite disks and found that variations of 
flcrit are < 0.5 AU. Residual heat from the moons’ accretion, ra¬ 
diogenic decay, and tidal heating might have provided additional 
heat sources, which we neglected. Hence, aciit = 4.5 + 0.5 AU 
must be considered a lower limit. 

If the GT actually took place and our conclusions about 
the pre-GT formation of the Galilean moons are correct, their 
early thermal evolution needs to be readdressed. Radiogenic de¬ 
cay in the rocky components and residual heat from accretion 
have been considered the main heat sources that determined 
post-accretion internal differentiation (Kirk & Stevenson 1987). 
The total heat flux of several TW in both Ganymede and Cal- 
listo (Mueller & McKinnon 1988) translates into a few tens of 
mW m“^ on the surface. Our results, however, suggest that sun¬ 
light might have been a signihcant external energy source. At 
1.5 AU from the Sun, illumination might have reached several 
W m“^, if the CSD was at least partly transparent in the optical. 
Near-surface temperatures might have been approximately 10 K 
higher for > 10^ yr (e.g. in Fig. 3a of Nagel et al. 2004). The 
GT might have significantly retarded the cooling of the Galilean 
satellites. 

Alternatively, Ganymede and Callisto might have become 
H 20 -rich after the GT, maybe through ablation of newly accreted 
planetesimals (Mosqueira et al. 2010), but then a mechanism is 
required that either prevented lo and Europa from accreting sig¬ 
nificant amounts of icy planetesimals or that triggered the loss 
of accreted ice. Tidal heating and the release of large amounts of 
kinetic energy from giant impacts might account for that. 

As H is an effective UV absorber (Glassgold et al. 2004) 
CSD gas might have shielded early Galilean atmospheres quite 
effectively. The net UV blocking effect depends on the disk scale 
height, the gas density prohle in Jupiter’s gap, and the residual 
gas flow (Fung et al. 2014). Ultraviolet photons might have been 
scattered deep into the gap, maybe even via a back-heating ef¬ 


fect from the dusty wall behind the gap (Hasegawa & Pudritz 
2010). This aspect of our theory needs deeper investigation and 
dedicated CSD simulations with Jupiter migrating to 1.5 AU. 


4. Conclusions 

In this Letter we show that the Grand Tack model for the mi¬ 
gration of Jupiter and Saturn, if valid, imposes important con¬ 
straints on the formation of their massive icy moons, Ganymede, 
Callisto, and Titan; (1) Ganymede and Callisto (probably also lo 
and Europa) formed prior to the GT (Sect. 2.1), (2) their forma¬ 
tion took place beyond 4.5 + 0.5 AU from the Sun (Sect. 2.1), (3) 
the Galilean moons would have lost any primordial atmospheres 
during the GT via non-thermal XUV-driven escape due to the ac¬ 
tive young Sun (Sect. 2.2), and (4) Titan’s thick N 2 atmosphere 
and constraints from moon migration in CPDs suggest that Titan 
formed after Saturn’s tack (Sect. 2.3). 

Detailed observations of the Galilean moons by ESA’s up¬ 
coming JUpiter ICy moons Explorer (JUICE), scheduled for 
launch in 2022 and arrival at Jupiter in 2030 (Grasset et al. 2013), 
could deliver fundamentally new insights into the migration his¬ 
tory of the giant planets. If Ganymede and Callisto formed prior 
to Jupiter’s Grand Tack, then JUICE might have the capabilities 
to detect features imprinted during the moons’ journey through 
the inner solar system. 
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